Strong spin relaxation length dependence on electric field gradients 



N£> 
O 
O 



ON 



c3 



13 
O 

o 



o 
o 

-I— > 
c3 



c 

o 
o 



D. Csontos 1 ' 2 ^ and S.E. Ulloa 2 
1 Institute of Fundamental Sciences, Massey University, 
Private Bag 11222, Palmerston North, New Zealand 
2 Nanoscale and Quantum Phenomena Institute, and Department of Physics and Astronomy, 
Ohio University, Athens, Ohio 45701-2979, USA 

We discuss the influence of electrical effects on spin transport, and in particular the propagation 
and relaxation of spin polarized electrons in the presence of inhomogeneous electric fields. We show 
that the spin relaxation length strongly depends on electric field gradients, and that significant 
suppression of electron spin polarization can occur as a result thereof. A discussion in terms of 
a drift-diffusion picture, and self-consistent numerical calculations based on a Boltzmann-Poisson 
approach shows that the spin relaxation length in fact can be of the order of the charge screening 
length. 

PACS numbers: 72.25.Dc, 72.25.Hg, 72.25.Rb, 72.25.Mk 

I. INTRODUCTION 

Electrical spin injection, transport and detection in diluted magnetic semiconductor (DMS)-nonmagnetic semicon- 
ductor (NMS) based systems^ are three crucial prerequisites for the successful realization of semiconductor spin- 
tronics. While the spin degree of freedom naturally is important, the charge degree of freedom plays a significant role 
in semiconductor spin transport as well, through electric field induced effects. Several experiments^^ have recently 
shown that the spin injection efficiency is strongly dependent on both applied and intrinsic electric fields, caused by, 
e.g., inhomogeneous doping. This was theoretically emphasized by Yu and Flatte(YF) 8 within a drift-diffusion model 
from which it was shown that the magnitude and sign of an electric field strongly influences the spin relaxation length 
according to 
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where E is a homogeneous electric field, and L s is an intrinsic (electric field-independent) spin diffusion length. 
According to eq. Q), the electric field enhances (suppresses) the spin relaxation length Lr}(Ljj) in the direction anti- 
parallel(parallel) to the direction of the electric field. However, the underlying assumption leading up to eq. (JTJ was 
local charge neutrality, and thus homogeneous electric fields. The question is how spin polarized electrons propagate 
in, e.g., inhomogeneously doped semiconductors where electric fields are inhomogeneous? 

In this paper, we use a similar approach as YF as well as self-consistent numerical calculations to study the effects 
of inhomogeneous electric fields. We find that, in an analogous drift-diffusion picture, a " quasi- local" spin relaxation 
length in the presence of inhomogeneous electric fields can be defined according to 
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from which a strong electric field gradient dependence becomes evident. However, in order to understand the full 
influence of inhomogeneous electric fields one needs to have access to the self-consistent charge density and electric field 
profile. We have performed such self-consistent numerical calculations using a Boltzmann-Poisson equation approach 
for a DMS-NMS inhomogeneously doped semiconductor from which we demonstrate that electric field gradients 
indeed can play a crucial role in the propagation of spin polarized electrons. Our findings may in particular highlight 
a possible difficulty in obtaining efficient spin polarized propagation at semiconductor interfaces between low-to-high 
doping concentration regions. 

In the following we will show a derivation of the drift-diffusion model leading to eq. (2) (Section 2), followed by 
our numerical model (Section 3), and a discussion of our results (Section 4). 
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FIG. 1: Schematics of the inhomogeneously doped DMS-NMS structure studied in our calculations. The DMS extends for 
x < x p , whereas the doping interface occurs at x = Xd- 



II. DRIFT-DIFFUSION DESCRIPTION OF SPIN TRANSPORT 

The current for spin-up and spin-down electrons can be written as 

JTU) = CT T(i) E + e£l T(i) Vn T(l) > ( 3 ) 

where -Dj(x) and c-m) are the diffusion constants and conductivities for the spin-up(down) species. Here, the con- 
ductivities refer to ct(n = e^tCD^tQ)' m terms of the mobilities "t(x) and densities 7if(n of spin-up(down) electrons. 
The continuity equations for the two species require that 

dn Ul) n Ul) | n l(T) , l v . r ,s 
oT = 1 H -V ■ Jt(i) > ( 4 ) 

dt T T1(1T) T lT(Ti) e 

where Tj-j, (t^) are the spin-flip scattering rates. In steady-state eqs. ij.'il 111 result in 

n entdi en iff) , 

W m) • E + o- TU) V • E + eD m) V 2 n m) = + _ (5) 

r U(it) T 1T(T1) 

For a NMS u j = /ij , and D j = £) j . Multiplying eq. JSJ with ctj, (o^ ) and substracting them from each other we arrive 
to the following expression for the spin density imbalance 5^ = n-f — rij 

V 2 % + ^V % + ^ u V.E-^0, (6) 

where we have used the Einstein relation u = eD/ksT, t s is the spin relaxation time r s _1 = tT, + r^j 1 (tjj = tjj), 
and L s = \JDr s is the intrinsic spin diffusion length. 

Equation 10 is a drift-diffusion equation for the spin density imbalance and contains terms depending on both 
the electric field, and the electric field gradient. We note that in fact eq. I© resembles the corresponding equations in 
Refs£, but with an additional term proportional to V ■ E. Similarly to Refs£ we can use the roots of the characteristic 
equation for eq. I© to define up- and down-stream spin diffusion lengths from eq. © leading up to eq. Notice 
that eq. © is only defined "locally", over a region where V • E can be considered constant, and where an average 
value of the electric field is used. Note also that the above equations are also valid for ferromagnetic semiconductors 
provided the mobility and diffusion constant are replaced with the ones corresponding to the minority-spin species^. 

A comparison between eq. © and Refs£ results in that eq. (0 above contains an additional term — eVS/fceT in 
the square root. Naturally, this term plays an important role in inhomogeneously doped semiconductors. We have 
performed self-consistent calculations and compared them with numerical fits using eq. (0) to verify its applicability. 



III. SELF-CONSISTENT NUMERICAL CALCULATIONS 



Our self-consistent approach is based on the solution of the Boltzmann and Poisson equations using numerical 
methods developed in Ref. 10 . Working in a one-dimensional model (choosing the x direction), the transport of 
spin-polarized electrons is described by two BTE equations according to 

eEdf ul) | ^ d/t(i) = -fr(i) ~ /"(j) /T(i)-/l(T) ^ 
m* dv dx T m T U(iT) 

where /t(j.) is the electron distribution for the spin- up (down) electrons, r m is the momentum relaxation time, and 
i/ T tl (1/Tj_t) is the scattering rate for spin-up (down) electrons. The first term on the right-hand side of eq. (JJJ 
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describes the relaxation of each nonequilibrium spin distribution to a local equilibrium (spin-dependent), normalized 
electron distribution function = n ](\:){ v )\J 2 ^k B T ex P ( — mv 2 /2kBT), where T = 300 K is the lattice temperature. 
The last term in eq. describes the relaxation of the spin polarization. From the distribution function /j-m we 
calculate the local spin density according to n^n\(x) = J fm)(r, v)dv, the total charge density, n = and the 

spin density imbalance, Sn = n-\ — n^. In order to take into account inhomogencous charge distributions and electric 
fields we solve eq. {7J) together with the Poisson equation V • E = e(No — n-f — ni)/(eeo) where e is the dielectric 
constant and No is the donor concentration profile, using finite difference and relaxation methods 10 . 

IV. RESULTS AND DISCUSSION 

We have performed self-consistent, numerical calculations for an inhomogeneously doped DMS-NMS structure 
schematically depicted in Fig.[l]using T = 300 K, m* = 0.067m , r m = 0.1 ps, r s = 0.5 ns, Ni = 10~ 21 m -3 , L = 0.2 
/im and Vb = —0.3 V. 

At the interface between the two doping regions with concentrations Ni and N 2 a built-in field is expected due to the 
diffusion of carriers from the region with higher to the region with lower doping concentration. This is illustrated in 
Fig. 0where we show the electric field profile, as well as the gradient of the electric field (inset) around the interfacial 
regions, for different doping concentrations N 2 = 0.05, 0.5, 1, 5, 20iVi. The sign and amplitude of the electric field 
gradients naturally depend strongly on the doping concentration N 2 on a length scale of the order of the screening 
length. 

While it is not possible to use eq. (0 directly for a quantitative analysis of the spin transport properties, one 
can make preliminary estimates of the order of the spin relaxation length around a specific region using linear 
approximations of the electric field. For example, using a linear fit of E for —0.1 < x < 0.1 fim for the N 2 — 20Aq 
curve yields WE ss 2.8 V//im 2 . Furthermore, using a mean value of E ps 2.9 kV/cm and L s ps 1.8 /im we obtain 
from eq. @ a spin relaxation length value L' D ps 0.16 /im along the (charge) direction of transport. Similarly, it 
is found that L'jj ps 0.06 /im. In comparison, an evaluation of eq. ifTjl. i.e. without the electric field gradient term, 
using the same value of the electric field yields ps 36 /im, and Ljj ps 0.09 /im. Thus, the spin relaxation length in 
the direction of the charge flow is predicted to be strongly influenced by the inhomogcneous electric field around the 
doping interface. 

In Fig. EI the numerically calculated spin density imbalance, S^, and total charge density, n, are shown for the 
N 2 = 20Ni structure considered in the above estimates. The total charge density (dashed line) increases monotonically 
as expected within a screening length from the low to the high doping concentration region. The spin density imbalance 
(solid line) follows a similar increase, but is, however, significantly suppressed such that only a maximum of 0.5A^2 
is reached around the interface region, beyond which 5-^ decreases exponentially. The observed suppression is a 
consequence of the inhomogeneous electric fields. 

A crude estimate of the validity of the values for L' D can be obtained by comparing our numerical calculations with 
S^i = n(x) exp[—(x — xa)/L' D ]. In Fig.0we show the resulting fits for L' D = 0.16, 0.3, 1 and 36 //m, the first and latter 
values corresponding to the previously estimated values for L' D and Lb obtained from the numerically calculated 
electric field profiles. We see that the higher values of Lo indeed yield that &n ps n. In contrast, the Ln values 0.16 
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FIG. 2: Electric field and electric field gradient (inset) profile for different doping concentrations N2 = 0.05, 0.5, 1, 5, 20iVi. 
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FIG. 3: Total charge density (dashed line) and spin density imbalance (solid line) profiles calculated for N2 = 2QN\. Additional 
curves for —0.1 < x < 0.1 /jm correspond to Sfi ~ n(x) exp[— (x + 0.1)/Ld] for different Ld as described in the text. Inset 
shows n and 5^ for a larger region. 



and 0.3 /im, which are of the order of the screening length, yield better fits to the numerically calculated data. 

Our results and analysis shows that inhomogeneous electric fields can strongly influence the spin relaxation length, 
and that significant spin polarization suppression can occur as a result thereof. We have also high-lighted a possible 
mechanism that may inhibit propagation of spin polarized electrons from low-to-high regions of doping concentrations, 
such as encountered in future semiconductor spintronic devices. 
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